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This mini-review is dedicated to the 85th birthday of Prof. L. V. Keldysh, from whom we 
have learned so much. In this paper we study the potential and electron density depth profiles 
in surface accumulation layers in crystals with a large and nonlinear dielectric response such as 
SrTiOa (STO) in the cases of planar, spherical and cylindrical geometries. The electron gas can 
be created by applying an induction Do to the STO surface. We describe the lattice dielectric 
response of STO using the Landau-Ginzburg free energy expansion and employ the Thomas-Fermi 
(TF) approximation for the electron gas. For the planar geometry we arrive at the electron density 
profile n(x) oc (x + d )~ 12 / 7 , where d oc We extend our results to overlapping electron gases 

in GTO/STO/GTO heterojunctions and electron gases created by spill-out from NSTO (heavily 
n-type doped STO) layers into STO. Generalization of our approach to a spherical donor cluster 
creating a big TF atom with electrons in STO brings us to the problem of supercharged nuclei. It is 
known that for an atom with nuclear charge Ze , where Z > 170, electrons collapse onto the nucleus 
resulting in a net charge Z n < Z. Here, instead of relativistic physics, the collapse is caused by 
the nonlinear dielectric response. Electrons collapse into the charged spherical donor cluster with 
radius R when its total charge number Z exceeds the critical value Z c ~ R/a, where a is the lattice 
constant. The net charge eZ n grows with Z until Z exceeds Z* ~ (R/a) 9 / 7 . After this point, the 
charge number of the compact core Z n remains ~ Z *, with the rest Z* electrons forming a sparse 
Thomas-Fermi atom with it. We extend our studies of collapse to the case of long cylindrical clusters 
as well. 


I. INTRODUCTION 

In recent years, there has been growing interest in 
the investigation of ABO 3 perovskite crystals, which are 
important for numerous technological applications and 
show intriguing magnetic, superconducting, and multi- 
ferroic propertied^. Special attention 2 3 J is paid to het¬ 
erostructures involving SrTi0 3 (STO) which is a semi¬ 
conductor with a band gap of ~ 3.2 eV® and a large 
dielectric constant k ranging from 2 • 10 4 at liquid he¬ 
lium temperatures to 350 at room temperature. As with 
conventional semiconductors, SrTiOa can be used as a 
building block for different types of devices, with reason¬ 
ably large mobilitj®^. 

Many devices are based on the accumulation layer 
of electrons near a heterojunction interface involving a 
moderately n-type doped STO. For example, one can 
get an electron accumulation layer on the STO side of 
the GTO/STO heterojunction induced by the electric 
field resulting from the “polar catastrophe” in GdTi0 3 
(GTO]P (see Fig. [l]). The role of GTO can also be played 
by perovskites LaAlOj 2 !!], NdA10 3 , LaVO^, PrA10 3 , 
NdGaOj^, LaGaOd^~d LaTiO^ 11 ! and others producing 
the polar catastrophe'. One can accumulate an electron 
gas using a field effeciP 2 ®^. In Refs. Il5l and ITfil the au¬ 
thors accumulated up to 10 14 cm -2 electrons on the sur¬ 
face of STO using ionic liquid gating. Inside bulk STO 
d-doping can be used to introduce two accumulation lay¬ 
ers of electrons 17 4 ^. It is natural to think that the depth 
profiles of the potential and electron density inside STO 
have a universal origin in all these devices. 

Interface properties determine characteristics of all 



FIG. 1. (Color online) Schematic energy diagram of an accu¬ 
mulation layer in a moderately n-doped STO. The band gap 
is E s . Electrons (blue region) are attracted by an external 
induction Do. The characteristic width of the electron gas is 
d. In the bulk of STO the Fermi level ef is near the bottom 
of the conduction band (plotted by the dashed line) 


these devices. Not surprisingly, the potential and elec¬ 
tron density depth profiles in such devices have at¬ 
tracted attention from the experimental 2 ^ 23 and theo¬ 
retical points of vievff^ 2 ^ 2 ^. For example, experimental 
data show that electrons are distributed in a layer of 
width ~ 5 — 10 nm near the LaA10 3 /SrTi0 3 interface. 
Theoretical works that attempt to explain such behavior 
are based on microscopic numerical calculations. 

The goal of this minireview based on Refs. [29] and 
mi is to create a simple, mostly phenomenological ap¬ 
proach for describing the electron density depth pro- 
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files in STO. To account for the nonlinear dielectric re¬ 
sponse in ST O we use the Landau-Ginzburg free en¬ 
ergy expansion? 1 Electrons are almost everywhere de¬ 
scribed in the Thomas-Fermi approximation^. Although 
we mostly concentrate on accumulation layers in STO, 
the developed approach is applicable to KTaOj® and 
CaTiOj^ serving as the host media for accumulation 
layers as well. 

Our main result is a new form for the potential and 
electron density depth profiles in various planar STO 
structures due to the nonlinear dielectric response. In 
particular, for an accumulation layer in STO created by 
an external interfacial induction Do, we find an electron 
concentration n(x) that depends on the distance from 
the surface x as n(x) oc (x + d)~ 12//7 , where the width 
d decreases with D 0 as d oc D^ 7 '' 5 . It is shown in Ref. 
1301 that these relations seem to agree with experimental 
dat ci 21 |22 l. Using this basic solution for a separate ac¬ 
cumulation layer we study more complicated problems 
where accumulation layers overlap, e.g., the structure of 
GTO/STO/GTO multi-heterojunctionP^SJ. We calcu¬ 
late how electron density profiles evolve as a function of 
the the distance between two heterojunctions. Another 
planar problem arises when the accumulation layer is cre¬ 
ated in STO which is a part of the modulation doped 
structure of NSTO/STCf 3 ®®. In this case electrons spill 
out from the heavily n-type doped STO layer (NSTO) 
and the induction Do is self-consistently provided by the 
depletion layer in NSTO. 

The experience with above planar problems allows us 
to study the electron gas in STO created by the exter¬ 
nal charges or electric field with spherical and cylindrical 
symmetry. Such devices cab be realized by doping the 
bulk STO by generating oxygen vacancies at high tem¬ 
peratures. The vacancies either form along a network 
of extended de fect^ or assemble together to lower the 
system’s energjSHHO producing large positively charged 
donor clusters. 

Another way to more controllably create such a cluster 
is to “draw” a disc of charge by the atomic force micro¬ 
scope (AFM) tip on the surface of LAO/STO str ucture 
with the subcritical thickness for LaA 103 (I.AO, 17 3 
The potential caused by such a positive disc in the bulk 
STO is similar to that of a charged sphere. 

Let us consider a spherical donor cluster with radius 
R and charge Ze. At relatively small Z , there are Z 
electrons located at distances from nb/Z to the Bohr ra¬ 
dius nb from the cluster, which form a Thomas-Fermi 
“atom’® with it. Here b = h 2 /m*e 2 , m* ss 1.5m is 
the effective electron mass in STC® with m being the 
free electron mass. Since k is large, the electrons are far 
away from the cluster and the whole “atom” is very big. 
As Z increases, the electron gas swells inward to hold 
more electrons. However, we find that when Z goes be¬ 
yond a certain value Z c ( nb/Z is still much larger than 
R at this moment), the physical picture is qualitatively 
altered. Surrounding electrons start to collapse into the 
cluster and the net cluster charge gets renormalized from 


Ze to Z n e with Z n <C Z at very large Z. 

The effect of charge renormalization is not ne W 50 * 51 l 
For a highly charged atomic nucleus with charge Ze , the 
vacuum is predicted to be unstable against creation of 
electron-positron pairs, resulting in a collapse of elec¬ 
trons onto the nucleus with positrons emitted^!. This 
instability happens when Z > Z c with Z c — 170 > 1/a, 
where a = e 2 /fic ~ 1/137 is the fine structure constant. 
When Z exceeds Z* ~ 1/a 3 / 2 ~ 137 3 / 2 , the net charge 
saturates at Z* (see Ref. ED- In the condensed mat¬ 
ter setting, there are similar phenomena in narrow-band 
gap semiconductors and Weyl semimetal^^ as well as 
graphene 52 . In all these cases, the collapse happens be¬ 
cause the energy dispersion of electrons is relativistic in 
the strong Coulomb field of a compact donor cluster play¬ 
ing the role of a nucleus. In our work, however, the col¬ 
lapse originates from the strong nonlinearity of dielectric 
constant in STO at small distances from the cluster. In 
the case of a spherical donor cluster, this nonlinearity 
leads to the change of the attractive potential near the 
cluster from being oc 1/r to oc 1/r 5 , resulting in the col¬ 
lapse of non-relativistic electrons to the cluster. 




FIG. 2. (Color online) The number of collapsed electrons 
S and the renormalized net charge Z n e as a function of the 
original charge Ze. S is shown by the thin solid line (red), 
Z n is denoted by the thick solid line (blue), and the dashed 
line (black) is a guide-to-eye where Z n = Z. Z c denotes 
the critical value where electrons begin to collapse and Z* is 
the saturation point where Z n stops growing, (a) Collapse 
of electrons and charge renormalization for highly charged 
nuclei. S oc Z J at Z c -C Z <S Z* (see Ref. [53). (b) Collapse 
of electrons and charge renormalization for spherical donor 
clusters in STO. S oc Z 9 / 2 at Z c C Z < Z*. 

The phenomena of electron collapse and charge renor¬ 
malization in both heavy nuclei and our work are pre- 
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sented in Fig. [2] In our case, the first electron collapses 
at Z ~ Z c ~ R/a , where a is the lattice constant, and 
at Z 3> Z* ~ (I?/a) 9 / 7 , the net charge of nucleus 
saturates as ~ Z*. 

We use the Thomas-Fermi approximation to show how 
the electron gas collapses into the cluster at Z Z c and 
find the corresponding electron density profile of the new 
two-scale Thomas-Fermi atom (see Fig. [9]). 

The remainder of this paper is organized as follows. 
In Sec. m we define the model based on the Landau- 
Ginzburg theory for calculating the lattice dielectric re¬ 
sponse and describe the parameters STO. In Sec. m we 
use the Thomas-Fermi approach to calculate the electron 
accumulation layer for a single GTO/STO heterojunc¬ 
tion induced by the polar catastrophe. In the Sec. [TV| 
we consider the case of multiple heterojunctions where 
accumulation layers from different interfaces overlap. In 
Sec. [V]we consider the case in which the electrons spill 
out from a heavily n-type doped STO material (NSTO) 
to the moderately n-type doped one for both a single 


junction and multi-junctions. In Secs. VI and |VII| we ex¬ 
tend our studies to the non-planar structure consisting of 
spherical or cylindrical donor clusters inside STO. In Sec. 
vraj we go beyond the zero temperature and investigate 


the thermal ionization of spherical donor clusters and its 
experimental implications. Sec. |IX| provides a summary 
and conclusion. 


II. THE MODEL 

Bulk STO typically is an n-type semiconductor with 
a concentration of donors N > 10 1 ' cm -3 . Let us dis¬ 
cuss the position of Fermi energy £p in such crystals. 
The electron spectrum near the bottom of the conduc¬ 
tion band is complicated'^, and in order to make the 
problem of an accumulation layer tractable analytically 
we assume that it is isotropic and non-degenerate with 
the effective mass m* ~ 1.5 where m is free elec¬ 
tron mass. Within the hydrogenic theory of shallow 
donors, the donor Bohr radius is equal to where 
b = h 2 /m*e 2 ~ 0.35 A, e is the electron charge, and k is 
dielectric constant of the material. At room temperature 
when k = 350, the Bohr radius nb = 123 A is so large 
that the Mott criterion for the metal-insulator transition 
in doped semiconductors N c b 3 = 0.02/k 3 leads to a very 
small critical concentration N c = 1-10 16 cm -3 . At helium 
temperatures k = 2 • 10 4 and N c = 6 ■ 10 10 cm -3 . Thus, 
at the experimentally relevant concentration of donors 
N > 10 17 cm -3 , we are dealing with a moderately doped 
semiconductor in which the Fermi energy lies in the con¬ 
duction band of STO. On the other hand, due to the 
relatively high effective mass the bulk Fermi energy Ep is 
smaller than the bending energy of the conduction band 
bottom near the interface (see Fig. [l]). For example, 
for N = 10 18 cm -3 , the low temperature Fermi energy 
calculated from the bottom of the conduction band is 
£p ~ 4 meV, which can be up to 100 times smaller than 


the bending energy of the conduction band bottom in an 
accumulation layer for GTO/STO. Therefore, we assume 
below that the Fermi energy coincides with the bottom 
of the conduction band. 

We are interested in the electron distribution near an 
interface of STO. We consider the case when the axis x is 
directed perpendicular to the interface (plane x = 0 ) and 
lies along the [100] axis of a cubic crystal of STO. (In fact, 
STO changes symmetry from cubic to tetragonal at T ~ 
110AT, but the distortion is smalP^and can be neglected). 
An external induction Dq applied from the left (see Fig. 
[l]) is directed along the x axis. In that case the problem 
is effectively one-dimensional. If the charge density is 
denoted by p(x), then the potential depth profile <p(x) in 
the system is determined by the equations: 


dD „ „ „ dv> „ , , 

— = 47 rp, D = E + 47 tP, — = — E , ( 1 ) 

ax ax 

where D(x), E(x), P(x) are the electric displacement 
field (the induction), the electric field and the electric po¬ 
larization in STO. Equations ([I]) should be solved with 
proper boundary conditions. For example, for a single ac¬ 
cumulation layer the boundary conditions are D( 0) = Dq 
and </?(oo) = 0 . 

To solve the system (JTJ) one needs to know two mate¬ 
rial relationships E(P) and p(ip). Let us start from the 
lattice dielectric response E(P). STO is well known as 
a quantum paraelectric, where the onset of ferroelectric 
order is suppressed by quantum fluctuation^. 

A powerful approach to describe the properties of 
ferroelectric-like materials is based on the Landau- 
Ginzburg theory. For a continuous second-order phase 
transition the Landau-Ginzburg expression of the free en¬ 
ergy density F is represented as a power series expansion 
with respect to the polarization P: 

F = F 0 +^P 2 + \a^P A - EP , (2) 

where Fq stands for the free energy density at P = 0 and 
r is the inverse susceptibility r = 4tt/(k — 1 ) ~ 4tt/k. 
In this work 0 < r -C 1, Pq = e/a 2 is the characteristic 
polarization and a ~ 3.9 A 5 ^is the lattice constant. The 
coefficient A describes the non-linear dielectric response. 
For all estimates below we use A = 0.8 following from 
Ref. [30. The last term of Eq. ([2]) is responsible for 
the interaction between the polarization and the electric 
field E. In general F depends on the components of the 
vector P, but in the chosen geometry the problem is one- 
dimensional, and all vectors are directed along the x axis. 
The crystal polarization P is determined by minimizing 
the free energy density F in the presence of the electric 
field E , SF/SP = 0. This condition relates E and P, 


_ 47t „ A o 


( 3 ) 
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We note that E <C 47 rP and thus D = E + 47 rP ~ 47 rP. 
The induction D c at which the transition from linear 
to nonlinear dielectric response occurs can be found by 
equating the first and second terms in the expression 0 : 


D c — Pq 



(4) 


If D <C D c the dielectric response of STO is linear and 
one can use the simplified expression for the electric field: 



(5) 


For D D c the dielectric response of STO is nonlinear 
and one must instead use the expression: 





( 11 ) 


where C x = (225/8)tt 2 ~ 278, C 2 = 1125tt/ 8 ~ 442. For 
a nonlinear dielectric response we obtain the equation for 
the potential: 


dx 


d \ 

dx ejb J 


1 / 3 ' 


2 3 / 2 1 
37t 2 6 4 / 3 


A 1 / 3 


;/ 6 2 \ 


2/3 


Po 




3/2 


) V/b) 


( 12 ) 

With the same boundary condition we get the solution: 


ip(x) = C 3 - - 


b\ 


8/7 


i J A 2 / 7 


8/7 


(13) 


E = 


A 

(4 


D 3 . 


( 6 ) 


n(x) = C 4 


1 

63 




(14) 


Next one should specify p(ip), which depends on the spe¬ 
cific device of interest. 


III. A SINGLE ACCUMULATION LAYER 


where C 3 = [5 6 3 6 tt 12 /(7 8 2 3 )] 1 / 7 ~ 5.8, C 4 = 

[5 9 3 2 tt 4 2 6 /7 12 ] 1/7 ~ 1.3. 

The characteristic length d can be obtained using the 
neutrality condition (see Eq. 0). For a linear dielectric 
response this gives: 


In a single heterojunction the external induction Dq 
caused by the “polar catastrophe” on the interface at¬ 
tracts electrons with a three-dimensional concentration 
n(x) inside the accumulation layer of STO (see Fig. 0. 
Our goal is to find the electron depth profile n(x) and its 
characteristic width d. 

Due to electric neutrality the number of accumulated 
electrons has to compensate the external field Do, i.e. 


OO 

4ne J n{x)dx = Do- 


(7) 


To take into account the electron screening of the external 
field we use the Thomas-Fermi approaclP^Slin which the 
electron concentration n(x) and self-consistent potential 
profile <p(x) are related as —e<p(x)+p(x) = £f = 0 , where 

/*(*) = (* r^^Wx)?/ 3 ( 8 ) 

is the chemical potential of the electron gas. Thus, one 
can obtain the solution of Eqs. 0 by replacing p(x) with 
en(x) and using relations 0, 0. For a linear dielectric 
response we obtain the equation for the potential: 


fjp_\ _ 2 3 / 2 1 1 /^\ 3/2 

dx 2 \e/bJ 37 r 2 b 2 k \e/b) 


(9) 


We use the boundary condition ip = 0 at x —> 00 and get 
the solutioJ 57 t 


ip(x) = CljK 2 
b 


x + d J 


4 


( 10 ) 


<i = c 56 (2) 2 , Ws(?(0 1/5 , (15) 

where C 5 = [7r 2 225/2] 1 / 5 ~ 4. For a nonlinear dielectric 
response: 


d = Cob 


(j) 


2/5 


e/a 


2 \ 7/5 1 


A 3 / 5 ’ 


(16) 


where Cq = (16/7)(5 2 3 2 7r 11 ) 1 / 5 ~ 84. The induction D c 
at which the transition from linear to nonlinear dielectric 


response occurs can be found from equating Eqs. (15) 
and (16). This gives 


D =^~J~ 

c Vac "' 1 ’ 


(17) 


where C 7 = (2 21 7r 9 /7 5 ) 1 / 6 ~ 12, consistent with Eq. 0. 
For STO, the critical field D c depends on temperature: 
D c ~ 0.1 e/a 2 for helium temperature and D c ~ 0.7e/a 2 
for room temperature. 

The three dimensional concentration profile n[x) for 
the nonlinear dielectric response Eq. (14) is the main 


result of this section. Note that n(x) has a very long 
tail with a weak 12/7 power law dependence, which may 
lead to some arbitrariness in measurements of the width 
of the electron gas. Indeed, only 39% of electrons are 
located within the distance 0 < x < d near the interface 
and 68 % of electrons are located within 0 < x < Ad. In 
the calculation above we used a continuous model. Actu¬ 
ally, along the [100] axis STO is composed of alternating 
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TABLE I. Percentage of electrons in the UO 2 layer M of STO 
for Do = 2ne/a 2 , corresponding to a total surface density of 
0.5e/a 2 


M 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

% 

27.9 

14.4 

9.0 

6.2 

4.6 

3.5 

2.8 

2.3 

1.9 

1.6 


TiC >2 and SrO layers. The conduction band of SrTiC >3 
corresponds to the bands composed of mainly 3 d orbitals 
of Ti. Integrating n(x ) over each lattice cell in Table [I] 
we get a percentage of electrons in each of the 10 first 
TiC >2 layers of STO for the case D = 2ire/a 2 . 


One can see from Eqs. (Ill and ( |14[ ) that the tails of 
the electron depth profiles n(x) at x d do not depend 
on Dq and behave like 


n{x) =C 2 t^k 3 ( - 
0 6 \x 


(18) 


and 


1 /t \ 12/7 

"<*>=<** U) 


A 3 / 7 


r 


(19) 


for linear and nonlinear dielectric responses, respectively. 
Even for Dq D c when the electron distribution n(x) 
at moderately large x is described by dependence (141, 


at very large distances the polarization becomes smaller 
and the linear dielectric response takes over so that the 
n(x) dependence switches from Eq. (14) to Eq. ©■ 
This happens at the distance 


10 = b (I) 


7/30 


31/10 © 


2/5 


.7/10 


( 20 ) 


(, Xq = 360 nm and 20 nm for helium and room temper¬ 
ature respectively). Thus, the tail of n(x) is universal. 
For small Dq < D c the tail has the form n(x) oc aY 6 . For 
Dq > D c it has the form n(x) oc x~ 12 '‘ for x < xq and 
n(x) oc x~ R for x > Xg- 

On the other hand, one has to remember that our the¬ 
ory is correct only when n{x) is larger than the concen¬ 
tration of donors in the bulk of the material. 

Let us verify whether the Thomas-Fermi approxima¬ 
tion is applicable at x ~ d, i.e., k F d 1. Here 
k F = (37T 2 ) 1 / 3 n(0) 1 / 3 is the wavevector of an electron 
at the Fermi level. For D 0 <C D c 


k F d = C 8 k 2 / 5 


while for Do S> D c 


k F d — Cg 


A 2 / 5 \aJ 


2/5 / n \ 1/5 

L>0 


e/a 2 

e/a 


\ 2/5 / „ /„2 \ 3 /5 


Do 


( 21 ) 


( 22 ) 


10 ~ 7 e/a 2 < Dq < 40 e/a 2 for room temperature. For 
lower temperatures this interval is even larger. Thus, 
the Thomas-Fermi (TF) approximation is applicable for 
practically all reasonable induction Do- 

Let us now consider how the applicability of the TF 
approximation holds with growing x. The TF p aram eter 
n x Dx can be estimated with the help of Eqs. (18 1 and 
(191 and we then get 


n x ' 3 x ~ 



X < Xo , 



X > Xo ■ 


(23) 


Thus, the TF approximation is valid in the whole range 
of nonlinear and linear dielectric response and fails only 
at very large distance from the interface x = Kb. 


IV. TWO OVERLAPPING ACCUMULATION 
LAYERS 


In the section m we have investigated a single accu¬ 
mulation layer induced in STO. In the GTO/STO/GTO 
structui'(P^“ an accumulation layer forms near each 
interface and these two layers overlap with each other. 
When the STO layer is thick, one can expect that the 
two accumulation layers overlap weakly by the vanishing 
tails and the final electron distribution can be described 
as a simple addition of two accumulation layers given 
by Eq. (14). (The induction Dq on each interface is 
47 re/ 2 a 2 > D c and we are interested only in the region 
where the distance to the interface is smaller than xo- 
So we use the nonlinear dielectric response here.) How¬ 
ever, as the STO layer gets thinner, the overlap becomes 
stronger. Due to the nonlinear physics here, the electron 
density profile changes substantially. Below we study the 
density profile of the electron gas in the GTO/STO/GTO 
structure where the width of the STO layer is 2 L (see Fig. 

§■ 


GTO 

STO GTO 

c 

1 2L 


FIG. 3. GTO/STO/GTO structure with wide enough GTO 
layers and an STO layer of width 2 L. The GTO layers gen¬ 
erate Do = ‘Ipe/a? on each interface. 


where Cg = ( 5 3 3 3 7 r 3 / 2 4 ) 1//5 ~ 6 , Cg = 4/7(157r 3 ) 3 / 5 ~ To make the mathematics more compact, in Secs. |IV| 
23. One can see that k F d > 1 in the range of 2 • and [Vj we employ the dimensionless notations, in which 
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a X 3/2 , £ > 0. (24) 


Eq. (12) is rewritten as 

_d_ / dxV /3 = 
dZ\dZ) 

Here £ = x/b is the scaled distance from one interface, 
X = <p/(e/b) is the potential in the units of e/6, and 
a = 2 3 ' 2 (a/6) 4 / 3 H 1 / 3 /37r 2 . Also, from Eq. S, we know 


i(x) = 


C 


10 


6 3 e/6 


3/2 


(25) 


where Cio = 2 3 / 2 /37r 2 « 0.1. The electron density n(x) 
can be scaled as n(£) = A 1 ^ 3 (a/b) 4 ^ 3 n(x)b 3 and we then 
get 

n(0 = ax 3/2 , (26) 

which is the dimensionless form of the TF approximation. 
The width of the STO layer is scaled as 2 L = 2L/b. 
According to Eq. (241, we can get 

\ 3/4 

5i , (27) 


where g i is a constant arising from the integration. By 
integrating Eq. (27), we get 


8 


rx( 0) 
>x(L) 


dx 


'§ olx 3/2 +9i) 


3/4 


= d£ = L. (28) 


In this structure, each GTO/STO heterojunction pro¬ 
vides a fixed number of electrons to the accumulation 
layer inside STO which is e/2 per unit cell with the in¬ 
terfacial induction D 0 = 4ne/2a 2 . Using Eq. © , this 
gives the value of dx/dt; on the £ = 0 and £ = 2 L in¬ 
terfaces, which is oc E oc Dq. Due to the symmetry, the 
electric field is zero in the middle of the STO layer which 
means dx/dl; = 0 at £ = L. We can choose a value for 
gi and calculate x(A) and y(0) using Eq . (|27| . Then we 
can put boundary values of % into Eq. (|28[) and get the 
corresponding value of L(gi). Reversing L(g i), we find 
the function gi(L). Therefore, at any given value of L 
we can pin down g\ and use Eq. (27) to numerically get 


the whole electron profile inside STO. 

To realize this, we need to try various values of g\ and 
tune accordingly until we find the L we want. In this 
process, it is necessary to know what values of g\ are 
physically possible. It is obvious that the extreme values 
of g\ appear at L —» 0 and L —» oo. For the former case, 
electrons are almost uniformly distributed over the thin 
STO layer with dx/d£ ~ 0 and the electron density is 
approximately 1/2 a 2 L everywhere which gives the value 
of x inside STO by Eq. (26) (see Fig. [4]). From Eq. (27), 
we know at dx/dt; = 0, a bigger value of x corresponds 
to a smaller value of g±. Thus, in this case we get the 
minimum value of g± as 

r , ~ i5/3 

8 A 1 / 3 (a/6) 4 / 3 6 3 (0.5/a 2 L) 


9 1 


5a 2 / 3 


(29) 


Since L can be arbitrarily small (the lattice constant a is 
regarded as infinitesimal), we get in fact g i mm = — oo. 

In the latter case where the STO layer is very thick, 
the two interfaces are quite independent and tails of ac¬ 
cumulation layers barely overlap. So, near the center of 
the STO layer where dx/dt; = 0, the electron density van¬ 
ishes. This gives the maximum value of g\ as g 1 max = 0. 
In this case, the electron distribution is close to the sim¬ 
ple addition of two accumulations layers described by Eq. 
(fill) (see Fig. El). 


So we get the domain of g± as (—oo, 0]. This means 
we can choose whatever value for gi and get physically 
meaningful results. By trying different values of g i, we 
can find the electron density profiles at certain L that 
we are interested in. In Fig. [4j we show our results for 
electron density profiles at 3 different values of L. We see 
the evolution from an almost constant n(x ) at L = 2a 
around which the quantum criticality is observecpS to 
the one reminiscent of two weakly overlapping tails of 
accumulation layers described by Eq. (14) at L = 8a. 


o 

IT 

£ 



x/L 


FIG. 4. (Color online) Electron distribution at different val¬ 
ues of L which is the half-width of the STO layer (a is the 
lattice constant). The thick solid lines are the electron con¬ 
centration n(x) in the unit of n(0) where n(0) is the electron 
concentration at x = 0 (red: L — 8a; black: L = 4a; blue: 
L = 2a). When L is not very large, the electron density profile 
is substantially different fro m th e addition of two accumula¬ 
tion layers described by Eq. (141. 


It is easy to see that n(x) is quite flat near the mid¬ 
dle. This is actually an effect of the nonlinear dielectric 
response as we show below. Indeed, due to the symme¬ 
try of the structure, for the potential x(£)> derivatives of 
odd order inclu ding dx/dl; are all vanishing in the mid¬ 
dle. Using Eq. (24), the second-order derivative d 2 x/dl- 2 
is found to be proportional to dx/dl^ and is then also zero 
at £ = L while the fourth-order is nonzero. This explains 
why the density changes so slowly near the middle of the 
STO layer. We emphasize again that this “flattening” 
effect originates from the nonlinear dielectric response. 
When the response is linear, the differential equation is 
described by Eq. ([9]), according to which the second- 
order derivative d 2 xjd £ 2 is nonzero even at dx/d£ = 0. 
The density change near the middle is then faster. 























V. SPILL-OUT OF ELECTRONS FROM 
HEAVILY DOPED n-TYPE STO (NSTO) INTO 
STO 


7 


In the NSTO/STO junctions P ^ H SI, the interface is 
formed between a moderately n-type doped STO with 
the Fermi level close to the conduction band bottom (be¬ 
low we simply denote it as STO) and a heavily n-type 
doped STO (NSTO), which has a much higher Fermi 
level. We start from considering a single junction be¬ 
tween thick layers of NSTO and STO. In this case as a 
result of the original difference between Fermi levels elec¬ 
trons spill out into STO to create a common Fermi level 
and the total number of spilled electrons depends only on 
the doping level inside NSTO. Below, we study the elec¬ 
tron distribution for this scenario (see Fig. [5]) assuming 
that the donor concentration in NSTO is so large that 
dielectric response is nonlinear. (The linear case of such 
spill-out problems was first addressed by Frenkel- 57 !) 

Inside NSTO, the charge concentration is [no — n{x)\e 
where no is the density of the positive background charge 
inside the doping layer and n{x) is the electron density at 
distance x from the interface. According to the Thomas- 
Ferrni approximation, we have 

1(f) = “X 3 ' 2 - fs°- < 30 > 


Here n 0 = A 1 ^ 3 (a/b) 4 ^ 3 n 0 b 3 is the scaled background 
charge concentration. From Eq. (30), it can be derived 
that 


dx 

dH 


7 ax 


5/2 _ 


4n 0 x + 92 


3/4 


£<0 


(31) 


X and dx/dl ; are continuous at £ = 0 which satisfy both 
Eqs. (31) and (341, we then have 


92 = 4n 0 x(0) = 


4n 0 Xc 

d 8 / 7 


Together with Eq. (32), we then have 

7/12 

>-) 


d = 


(35) 


(36) 


which gives the expression for the decay length. 

Now we get the electron distribution inside STO. Using 
the differential equation (31) and the boundary values of 


X and dx/d^ at £ = 0, the electron density profile inside 
NSTO can also be obtained numerically. A schematic 
plot of the electron distribution is presented in Fig. [5] 
This is the universal curve independent of the value of 
n 0 - 



x/d 

FIG. 5. The electron distribution in a single NSTO/STO 
junction. The distance x is in units of the decay length d = bd. 
The electron density n(x) is in units of the donor concentra¬ 
tion no- This figure is universal for all values of no¬ 


where g 2 is a constant that we determine from the bound¬ 
ary conditions. Since at £ —> — oo, n(£) = no, we have 
dx/dt, • = 0 and x = ino/a) 2 ^. Therefore, 


12 n. 


5/3 


92 = 


(32) 


5 a 2 / 3 

Now let us switch to the STO side. By rewriting Eq. 


(|13|) in the dimensionless form, we have in STO 

£>0 


which gives 


X = 


dx 


Xc 


(£ + d ) 8 / 7 


(33) 




\ 3/4 


(34) 


Here d = d/5 is the scaled decay length and Xc = 
(2 3 5 3 /7 4 a 3 )“ ‘. Since the interfacial field Dp is now no 


longer a fixed value, we cannot use Eq. (16) to get the 


decay length. Instead, using the boundary condition that 


Let us now dwell on the case of a heavily doped 
NSTO layer with the width 2t embedded in STO 
(STO/NSTO/STO structure). Experimentally used 6 
layers have sharp boundaries 58 with well defined volume 
concentration no. One can imagine two limiting cases. 
When NSTO layer is so wide that its half-width t is much 
larger than the width of the electron gas in STO, one can 
use the above single-junction theory for the two separate 
STO/NSTO and NSTO/STO junctions. 

In the opposite case when t is much smaller than the 
width d of the electron gas on each side, most of elec¬ 
trons are located outside of such d-layer. As a result 
the effective two-dimensional charge density 2 tno cre¬ 
ates two accumulation electron gases on both sides from 
the d-layer. The external electric field for each gas is 
D 0 = 47rnot. The width of the gas is determined by Eq. 
(16) if Dq > D c which is given by Eq. 


Below we discuss how such peculiar distribution of elec¬ 
trons may affect the low temperature transport proper¬ 
ties of the structure. Namely, for a fixed concentration no 
we find the ratio of the mobility fig for a d-doped struc¬ 
ture with t d to the bulk mobility /tf, for a structure 
with t d. 


















We focus on low temperatures where electrons are 
mostly scattered by donors and the phonon scattering 
is negligible. The scattering crossection of an electron 
on an ionized donor is E ~ a 2 (a is the lattice constant) 
and does not depend on the electron energy. Indeed, the 
Coulomb impact parameter for an electron with Fermi 
energy e F = (3n 2 n) 2 ^ 3 h 2 /2m is I = e 2 /k£ F ■ For con¬ 
centration n larger than 10 1 ' cm -3 , I -C a. As a re¬ 
sult the cross-section saturates at the neutral-impurity 
cross-section E ~ a 2 . This unique feature of materi¬ 
als with high dielectric constant!^ simplifies the calcu¬ 
lation of the mobility. In the case of (5-doping, most of 
electrons live outside of the layer and are scattered by 
donors only when they move through this layer. This 
leads to a reduced effective concentration of scattering 
centers not/d which is much smaller than the concentra¬ 
tion no in the t d case. At t <C d the Fermi velocity 
of the electron gas is vs — {h/m)(jiot/d) 1 ^ 3 oc t 4/5 

since 

d oc t -7 / 5 for the nonlinear dielectric response, while for 
t^$>d the velocity vs does not depend on t. In the result, 
Us/nb = ( d/t)(vb/vs ) oc f -16 / 5 . Thus, the mobility fj,$ 
decreases with t at t <C d and at t ~ d saturates as the 
bulk value This effect probably was observed in Ref. 

1S01 

Let us now address more complicated periodic struc¬ 
tures formed by NSTO and STCfSHMI ( see Fig. |@j). Inside 


STO 

NSTO 

STO 

NSTO 


-2t 


0 


2L 


FIG. 6. Periodic NSTO/STO structure. The width of the 
NSTO layer is 2 1 and the width of the STO layer is 2 L. 


each NSTO layer, the potential should satisfy Eq. (31) 


while in STO, it obeys Eq. ( |27| . Similar to Eq. ([28]), we 
get 


rx(-t) 


d'X 


'x(o) (fay 5 / 2 - 4n 0 X + 32) 


3/4 


= J = t. (37) 


where t = t/b and 21 is the dimensionless width of the 
NSTO layer. Similarly to what we have done above, we 
can choose certain values for 32 and x(0) = Xo w here x(0) 
is the scaled potential on the interface. Then we can cal¬ 
culate corresponding t and L and the electron density 
profile. Thus, in reverse, at any given t and L we can 
find the electron distribution. Again, we need to find 
the physical range of g 2 and vn, w hich can be done sim¬ 
ilarly to what we did in Sec. IV Given values of t and 
L, we try different <72 and xo from their domains until 
we find the t and L we want. In this way, we get the 
electron distribution for 3 different values of L in Fig. [ 7 ] 


(for simplicity, we choose t = L here). We see that as 
L = bL increases, the periodic electron density profile 
n(x) evolves from being relatively constant to strongly 
oscillating. Again, one can see a “flattening” effect near 
the middle of each layer, either STO or NSTO. Like what 
we did in Sec. |IV[ one can verify that all derivatives of \ 
vanish in the middle of both layers until the fourth-order 
one. 



FIG. 7. (Color online) Electron distribution from —t to L 
plotted at t = L for different values of L. The thick solid 
lines are the electron density n(x ) in the unit of no where no 
is the donor concentration in NSTO (red: L = 10a; black: 
L = 6a; blue: L = 2a). This graph is plotted at no = 0.5/a 3 
where a is the lattice constant. As the layers get thinner, the 
electron distribution gets more uniform, deviating from the 
one shown in Fig. [5] 


Let us also dwell on a recently discovered new type of 
STO-based heterojunction formed by STO and NdTiOa 
(NTO). In this structure, the ceiling of the valence band 
created by the Mott gap in NTO is above the bottom of 
the conduction band of STCjhlFAf So in addition to the 
“polar discontinuity” 5 &2MI, the broken-gap band align¬ 
ment further brings electrons into STO by the electron 
’’spill-out” and the electron distribution becomes more 
complicated. One can basically use the methods we em¬ 
ployed in this section to solve the problem. The differ¬ 
ence due to the presence of the “polar discontinuity” is 
the abrupt jump of the induction field on the interface 
which gains an additional value of 47re/2a 2 when going 
from NTO to STO. One should note that the dielectric 
response inside NTO is linear and the potential obeys 
Eq. ([9]). Therefore one gets an expression of the poten¬ 
tial derivative dx/d£ different from Eq. (31). Also one 
should know the density of states below the Mott gap in 
NTO in order to solve for the electron distribution. The 
rest of the procedure is quite similar and one can then 
obtain the electron density profile in this structure. 


VI. SPHERICAL DONOR CLUSTERS 


So far, we have considered only the planar structures 
based on STO. In Secs. [Vl]and[VnJ we extend our stud¬ 
ies to other non-planar geometries (spherical and cylin¬ 
drical). To make the physics clearer, we return to the 
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dimensional expressions for physical quantities from now 
on as in Sec. cm 


A. Renormalization of Charge 

Consider a large spherical donor cluster of the radius R 
and the total positive charge Ze such that a <C R < nb/Z 
(for example, R can be 3 nm and Z ~ 60). If the di¬ 
electric response is linear, the electrons are mainly lo¬ 
cated at distances between rq = refe/Z and r a = nb from 
the cluster®^. For a very large re, these radii are huge 
([ta = 700 nm in STO at liquid helium temperature) and 
the electrons are far away from the cluster. However, at 
small distances, the dielectric response is nonlinear and 
changes the potential form. If the potential energy out¬ 
weighs the kinetic energy, electrons are attracted to the 
cluster and renormalize the net charge. To see when this 
happens, we look at the specific form of electric potential 
in this situation. We can calculate the potential from the 
differential equations (38a) and (38b) applicable to the 
spherical structure: 


d_ 

dr 

d_ 

dr 


dip \ 1 / 3 A 1 / 3 e . . . . . . 

= — 2/3 H r ) ~ n o],r<R (38a) 


^v /3 = 

dr ) 


P 
M) 

A 1 / 3 e 


P: 


,2/3 


■ 0 ), 


r> R (38b) 


o 


where r is the radius from the cluster center, n(r) is the 
electron density at radius r and no is the donor concen¬ 
tration inside the cluster. However, due to the simple 
charge distribution here, we can get the potential in an 
easier way. At r > R, the sphere looks like a point charge 
and D(r) = Ze/r 2 . Using this together with Eqs. |l]) 
and one can calculate the electric field and get the 
electric potential p(r) as: 


p(r) = 


A 

p2 

M) 


Ze 

47T 


1 

5 r 5 ’ 


R < r <C rq (39) 


no bound states of electron in the cluster. However, when 
Z is big enough so that Z > Z c , the electron can have 
negative total energy at r < R and will collapse into the 
cluster. Using Eqs. (39) and (|40|, we find 


a a 


(41) 


As Z continues increasing, more and more electrons get 
inside the cluster filling it from the center where the po¬ 
tential energy is lowest (see Fig. [8]). The single-electron 
picture no longer applies. Instead, we use the Thomas- 
Fermi approximation 49 with the electrochemical poten¬ 


tial p = 0, which gives Eq. (25). (We continue to assume 
here that the bulk STO is a moderately doped semicon¬ 
ductor.) 

When the number of collapsed electrons S is small, 
their influence on the electric potential is weak. One 
can still use Eqs. (39) and (40) for tp(r) and get the 


corresponding expression of n(r) 7 At r > R : since p(r) is 
oc 1/r 5 , we get n{r) oc I/ 7 ' 15 / 2 . In this way, we calculate 
S as 


S = 


where 


p°o ( Z \ 7 2 

J n(r)4:TTr 2 dr = 0.5 Z f — j oc Z 9 / 2 , 


Z* = 


Aitib/Aaf^R 


9/7 


The net charge number of the cluster is 


z n = z-s = z 


1 - 0.5 ( 

z* 


7/2' 


(42) 


(43) 


(44) 


One can see, when Z c <C Z « Z*, one gets S<Z and 
Z n Z, meaning the charge renormalization is weak. 
However, at Z ~ Z*, according to Eqs. (42) and (441, 
we get Z„ ~ S ~ Z*. The potential contributed by 
electrons is no longer perturbative. This brings us to the 
new regime of strong renormalization of charge. 


with ip(r = 00 ) defined as zero. Inside the cluster at 
r < R, since the charge is uniformly distributed over the 
sphere, the total positive charge enclosed in the sphere 
of radius x is equal to Zer 3 /R 3 , so D(r) = Zer/R 3 . One 
then gets the corresponding potential p(r): 


ip{r) 


A_ / Ze\ 3 f9_J_ 

P^\4^J \20 i? 3 


4 I? 9 J 


, 0 < r < R 


(40) 


using the boundary condition ip(r = R~) = ip{r = R + ). 
A schematic graph of the potential energy U(r) = —ep{r) 
is shown in Fig. [8] by the thick solid line (blue). 

The Hamiltonian for a single electron is H = p 2 /2m* — 
eip(r), where p is the momentum of the electron and m* is 
the effective electron mass in STC 30 . If we approximately 
set p ~ h/2r, we get a positive total energy of the electron 
everywhere when Z is very small. This means there are 


We show that at Z > Z* the net charge Z„e satu¬ 
rates at the level of Z*e. Indeed, when Z grows beyond 
Z*, Z„ can not go down and therefore can’t be much 
smaller than Z*. At the same time it can not continue 
going up, otherwise as follows from Eqs. (39) and (42) 
with Z replaced by Z n Z*, the total electron charge 
surrounding the charge Z„e at r > R would become 
Se ~ Z„e(Z n /Z*) 7 / 2 /g> Z n e leading to a negative charge 
seen from infinity. Thus, at Z /§> Z*, the net charge 
Z n saturates at the universal value of the order of Z*. 
This result is qualitatively similar to the one obtained 
for heavy nuclei and donor clusters in Weyl semimetals 
and narrow-band gap semiconductors in Ref. [5l] 

In the following subsection, we show how the renormal¬ 
ization of charge at Z Z* is realized through certain 
distribution of electrons, in which a structure of “double 
layer” (see Fig. |8| plays an important role. 
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to the surface and the characteristic decay length d <C R 
is given by Eq. (361 which in the dimensional form is 


d = 


C 


n 


A 1 / 4 


1/4 a 
a) (n 0 a 3 ) 7 / 12 ’ 


(45) 


where Cn « 2. By expressing no in terms of Z and R, 
we get d/R oc ( Z*/Z ) V12 <1 at Z > Z*. 

Correspondingly, the radial electron concentration at 
x > R is given by 


FIG. 8. (Color online) Potential energy of electrons U(r) = 
—eip(r) caused by the spherical donor cluster of radius R as 
a function of distance r from the cluster center, ip o is defined 
as n(<po) = no, where no = 2>Z/A-kR 3 , n(r) is a function of 
p(r) given by Eq. ( |25| . The thick solid line (blue) represents 
the potential profile of a cluster of charge Z < Z* which is 
in the regime of weak charge renormalization. The thin solid 
line (red) represents the potential of a cluster at Z > Z* 
in the strong renormalization regime, where the two verti¬ 
cal dotted lines show edges of the “double-layer” structure of 
width ~ U. The horizontal dashed line (black) indicates 
the position of the chemical potential p = 0. r* is the exter¬ 
nal radius of the collapsed electron gas where Thomas-Fermi 
approach fails. 


B. Radial Distribution of Electrons 


At Z Z*, the charge renormalization is strong and 
the most of the sphere of radius R is completely neutral¬ 
ized by electrons. In the neutral center of the sphere, 
the electron density n(r ) = no, where no = 8 Z/AttR 3 is 
the density of the positive charge inside the cluster. The 
corresponding “internal” electric potential <Pi n (r) = po 
where <po is given by n(<po) = no using Eq. (25). v 5 m( r ) 
is then oc (noa 3 ) 2 / 3 oc [Z/(R/a) 3 ] 2 ^ 3 . Outside the clus¬ 
ter, when the charge is renormalized to Z n , one gets a 
potential p 0 ut(x) similar to Eq. ( [39] ) with Z replaced by 
Z n . Since Z n is ~ Z* where Z* is given by Eq. (43), 
we get (fiout{r) oc (R/a)~ 8 ^ 7 at a distance r of the order 
R. Thus, close to the cluster surface, the ratio of the 
outside potential Poutix) to the inside potential Pi n (r) is 
~ ( R/a) 6 / 7 /Z 2 / 3 < 1 since Z > Z* ~ (. R/af I 7 . This 
indicates a sharp potential drop across the sphere surface. 

At 0 < R — r <C R, there’s a thin layer of uncom¬ 
pensated positive charges. At0<r — R R, a higher 
potential than farther away means a larger electron con¬ 
centration that forms a negative layer close to the sur¬ 
face. This “double-layer” structure resembles a capacitor 
which quickly brings the potential down across the sur¬ 
face as shown in Fi g. [51 An analogous structure also 
exists in heavy nucleP2E3J with charge Z 1/a 3 / 2 . 

To make the analysis more quantitative, one needs to 
know the specific potential profile in this region. Near the 
cluster surface, we can approximately use a plane solution 
of (M r), i.e ., ign ore the 2 jr term on the left side of Eqs. 
Bal and (38b). This kind of solution for r > R is given 


by Eq. (13) with x = r — R -C R which is the distance 


n(r)r 2 = 


C*12 1 

(43/753 

C \2 1 

(43/753 


b\ 


b\ 


24/7 

( a N 

12/7 

( 

\x + d j 

) 

24/7 

I 

( a N 

12/7 

1 

\x + d , 

) 


(46) 


R 


where r « R, C 12 ~ 1. 

Since the “double-layer” structure resembles a plane 
capacitor, near the surface, the potential drops practi¬ 
cally linearly with the radius. Using Eq. (13), one can 
get p(r) ~ [1 — 8x/7d]p(R) atO <x = r — R<^Ld, 
which gives the electric field 8p(R)/7d inside the “dou¬ 
ble layer”. At r < R, this electric field persists and gives 
ip(r) « [1 + 8(— x)/7d]ip(R) at 0 < —x = R — r <C d. 
As r further decreases, the positive layer ends and the 
potential crosses over to the constant value po given by 
n(ipo) = no using Eq. (25). 

According to Eq. (13), when x = r — R is comparable 


to R and the plane approximation is about to lose its va¬ 
lidity, p(r) is oc ( R/a )~ 8 / 7 . It is weak enough to match 
the low electric potential yw(r) a ( R/a )~®/ 7 caused by 
the renormalized charge Z n ~ Z* at r ~ R. The plane so¬ 
lution then crosses over to the potential ip 0 ut(r) oc Z * 3 /r 5 
which is the asymptotic form at large distances. 

A schematic plot of the potential energy U(r) = 
—e<p(r) as a function of radius r is shown in Fig. [|by 
the thin solid line (red). The corresponding radial distri¬ 
bution of electrons is shown in Fig. [9] by the thick solid 
line (red). 


So far, we have got a 1/r 5 potential p(r) and 1/r 11 / 2 
radial electron concentration n(r)r 2 at r R in both 
weak and strong charge renormalization cases. However, 
as the electron density decreases to certain extent so that 
the Fermi wavelength A is comparable to the radius r, 
the gas is no longer degenerate and the Thomas-Fermi 
approach fails. Since A ~ n(r) _1 - /3 , we get this radius 
r* ~ Z*a at Z 3> Z*. One should then return to the 
Schrodinger equation used for a single electron. Since the 
uncertainty principle estimates that the kinetic energy 
decays as 1/r 2 while the potential energy is oc — Z* 3 jr 5 , 
the potential energy is smaller than the kinetic energy 
in magnitude at r > r*, which means electrons can not 
stay at radii larger than r*. One can also find that us¬ 
ing the Thomas-Fermi solution ip out (r), the total electron 
number calculated at r > r* is ~ 1, which again indicates 
there is no electron at r > r* considering the discreteness 
of electron charge. As a result, the 1/r 11 / 2 tail of radial 
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FIG. 9. (Color online) Radial electron concentration n(r)r 2 as 
a function of radius r for the two-scale Thomas-Fermi atom 
formed by the super-critical cluster of donors (solid lines). 
The thick solid line (red) represents the inner collapsed elec¬ 
trons at r < r c where the dielectric response is nonlinear. 
The thin solid line (blue) shows the electrons belonging to 
the outer shell which form the standard Thomas-Fermi atom 
with the renormalized nucleus of charge Z* at r > r c , where 
the dielectric response is linear. This electron gas ends at the 
Bohr radius rA = while most of them are at the radius 
r m = nb/Z* 1 ^ 3 . For contrast, the dashed line (black) den otes 
the electrons forming a conventional Thomas-Fermi atonP 3 ® 
with a nucleus of charge Z when Po is infinity and there’s no 
range with nonlinear dielectric response. The reduction of 
electron density in the outer shell of electrons due to the col¬ 
lapse is substantial. The reason this is not immediately seen 
from the difference of height between the dashed line (black) 
and the thin solid line (blue) is that we use a logarithmic 
scale here, no is defined as 3Z/4nR 3 . This graph is plotted 
at b = 0.35 A, a = 3.9 A, A = 0.9, R = 4.4 a, k = 20000, 
no = 0.8/a 3 . 


electron concentration will not continue to infinity but 
stop at radius r*. This is a semi-classical result. Quan¬ 
tum mechanical analysis shows that the electron density 
does not go to zero right at r* but decays exponentially 
after this point. Since this decay is fast and brings very 
small corrections to the edge of the inner electron gas, 
we do not consider it here. 

At k = oo, the rest of the electrons are at the infinity 
so that we are dealing with a positive ion with charge Z*. 
At finite but very large k, at certain distance from the 
cluster, the field is so small that P > yj 4-7T / kAPq is no 
longer satisfied and the linear dielectric response is recov¬ 
ered. Things then become quite familiar. Electrons are 
mainly located between n = nb/Z* and = k b with 
the majority at radius r m = nb/Z * 1 ^ 3 as given by the 
Thomas-Fermi modeP^. Although quantum mechanics 
gives a nonzero electron density at r < rq, the num¬ 
ber of total electrons within this radius is only ~ 1 and 
can be ignored. So approximately, when rq r*, i.e., 
k (Z*) 2 ~ (P/a) 18 ' 7 , there’s a spatial separation be¬ 
tween inner collapsed electrons and outer ones that form 
the usual Thomas-Fermi atom with the renormalized nu¬ 
cleus. When k is not so big, such separation is absent, 
which actually happens more often in real situations. The 
inner tail then connects to the outer electrons with the 
Thomas-Fermi approach valid all the way and the dielec¬ 


tric response becomes linear at r = r c oc atA^Z* 1 ^ 2 . 
One should note, as long as k is large enough to sat¬ 
isfy r m R which gives k ( R/a ) 10 ' 7 , the majority 
of the outer electrons located at r m do not intrude into 
the cluster or the highly screening double-layer structure 
near the cluster surface. The charge renormalization pro¬ 
cess remains undisturbed and the total net charge seen 
by outer electrons is still Z* . The corresponding radial 
electron concentration n(r)r 2 is shown in Fig. [9] At 
k < (R/a) 10 / 7 , in most of the space the dielectric re¬ 
sponse is linear. In that case, almost all electrons re¬ 
side in the cluster with only some spill-out near the sur¬ 
face. The positive and negative charges are uniformly 
distributed inside the cluster as described by the Thomp¬ 
son “jelly” model. 


VII. CYLINDRICAL DONOR CLUSTERS 


In some cases, the donor clusters are more like long 
cylinders than spheres. Then, a cluster is described by 
the linear charge density i]e while its radius is still de¬ 
noted as R. We use a cylindrical coordinate system with 
the z axis along the axis of the cylinder cluster and r 
as the distance from the axis. We show that when the 
charge density rje is larger than certain value r] c e, elec¬ 
trons begin to collapse into the cluster and the charge 
density is weakly renormalized. When 77 exceeds another 
value 77* 77c, the renormalization becomes so strong 

that the net density rj n remains ~ rf regardless of the 
original density 77. Our problem is similar to that of the 
charged vacuum condensate near superconducting cosmic 
stringd^, and is also reminis cent of the Onsager-Manning 
condensation in salty watei^® 7 !. 

Renormalization of Linear Charge Density. — For a 
uniformly charged cylindrical cluster with a linear charge 
density r)e, similar to what we did in Sec. m we 
get D(r) = 2rj(r)e/r , where r](r) is the total linear 
charge density enclosed in the cylinder of radius r and 
r)(r) = rjr 2 / R 2 at r < R and 77(r) = 77 at r > R. We then 
can calculate the electric field using Eqs. 0 and 0 and 
get the electric potential <p(r) as: 


<p(r) 

<P(r) 


A /pe\ 3 / 3 1 
P 2 \2n) {Ir 2 
A / r]e\ 3 1 
P 2 V2W 2H’ 


1 r 4 \ 
4 Rp) 


, 0 < r < R 


( 47 a) 


R<r ( 47 b) 


with ip(r = 00) chosen to be 0 . The corresponding po¬ 
tential energy U(r) = —e<p(r) is shown in Fig. 11 by the 
thick solid line (blue). Using the Schrodinger equation 
and setting the momentum p — h/2r, we find that the 
tightly bound states of electrons, in which electrons are 
strongly confined within the cluster (at r < i?), exist only 
when 77 > 77 c where 


7?c « 27T 



( 48 ) 
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which, contrary to Z c obtained in the spherical case, does 
not depend on R. Electrons begin to collapse into the 
cluster at r] > rj c and in the beginning they are located 
near the axis where the potential energy is lowest (see 
Fig. 11). With increasing 77 , the electron density grows 


and one can adopt the Thomas-Fermi description. Using 
Eq. (25) and (47b), one gets the electron density n(r) oc 
1 /r 3 at r > R and the total number of collapsed electron 
per unit length is 


/*oo / \ 7/2 

0 = J n{r)2irrdr = 0.5r] J oc rf^ 2 , 


(49) 


where 


V = 


21,1 -L) 


1/3' 


9/7 


d\ 2/7 

f ' 


(50) 


of the characteristic decay length d is also the same as 
in Eq. ( |45| . When x = r — R is comparable to R, the 
plane solution crosses o ver to the fast decaying potential 
oc 1/r 2 as given by Eq. (47b) with 77 replaced by ?/„ ~ rf. 
A schematic plot of the potential energy U(r ) = —&p(r) 
is shown in Fig. m 



The net charge density r] n e is then renormalized to 

7/2" 


Vn = V — 0 = V 


1 - 0.5 ( -9- 

rf 


(51) 


At 77 <C 77 *, the renormalization of charge density is 
weak and r] n grows with 77 . At 77 > 77 *, the renormaliza¬ 
tion effect becomes strong. Most of the cluster is then 
neutralized by electrons and the final net density rj n is 
much smaller than 77 . Following the logics similar to that 
in the spherical case, and by using Eq. (49), one can 
show that rj n reaches a saturation value of 77 * at 77 77 *. 

The dependence of r] n on 77 is shown in Fig. [To] 



n‘(R) 


FIG. 11. (Color online) Potential energy of electrons U(r) = 
—eip(r) caused by a cylindrical donor cluster of radius R as 
a function of the distance r from the cluster center. <po is 
defined as n(ipo) = no, where no = r//nR 2 , n(r) is a func¬ 
tion of <p(r) given by Eq. ( |25| ) . The thick solid line (blue) 
represents the potential profile of a cluster of charge density 
77 < 77 * which is in the regime of weak renormalization of 
charge. The thin solid line (red) represents the potential of 
a cluster with 77 3 > 77* which is in the strong renormaliza¬ 
tion regime. The two vertical dotted lines show edges of the 
“double-layer” structure of width ~ d <C R. The horizon¬ 
tal dashed line (black) indicates the position of the chemical 
potential 77 = 0 . 

This potential produces a universal tail of electron den¬ 
sity n(r) ~ 1/r 3 . The corresponding radial electron con¬ 
centration n(r)r is ~ 1/r 2 . Since the Fermi wavelength 
A ~ 7 t(r) -1 / 3 , we get A ~ r, i.e., the Thomas-Fermi ap¬ 
proach is only marginally valid. The collapsed electrons 
extend until the linear dielectric response is recovered 
and then connect to the outer electrons. 


VIII. FINITE-TEMPERATURE EFFECT IN 
SPHERICAL DONOR CLUSTERS AND ITS 
EXPERIMENTAL IMPLICATIONS 


FIG. 10. (Color online) Number of collapsed electrons per 
unit length 9 and renormalized net linear charge density r/ n 
as a function of the cluster linear charge density rj. The thick 
solid line (blue) shows 7771 ( 77 ). The thin solid line (red) rep¬ 
resents #( 77 ). The dashed line (black) is a guide-to-eye with 
t ?71 = 77 . 9(rj) oc 77 ®' 2 at 77 c -C 77 < 77 *. 


Radial Distribution of Electrons .— At 77 77 *, there 

are lots of collapsed electrons inside the cluster where 
n(r) = n 0 = r]/nR 2 and the potential energy is low. 
Again, there’s a “double-layer” structure on the surface 
that provides steep growth of potential energy with r at 
r = R. Close to the cylinder surface atO <x = r — f?<C 
R, as for the sphere, we can approximately use a plane 
solution of tp(r) as given by Eq. (13). The expression 


So far, we dealt with zero temperature. At a finite 
temperature T, the neutral cluster atom can get ionized 
due to the entropy gain of ionized electrons. The donor 
cluster atom becomes a positive ion with charge Zi(T)e. 
Recall that the TF approach is valid at all distances until 
rA = nb. The first ionization energy of the cluster atom 
is then negligible (= e 2 / n 2 b ~ 10 ~ 7 eV). So the cluster 
atoms are always partially ionized. Our goal below is to 
find this ionization charge for the spherical donor clus¬ 
ters. Similar analyses for the cylindrical clusters can be 
found in Ref. [25] and are not repeated here. 

We assume that we have a small but finite three- 
dimensional concentration N of spherical clusters and 
the charge Zi(T) < Z*, i. e., the outer electron shell 
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is still incompletely ionized. Such a cluster can bind 
electrons with an ionization energy Zi(T) 2 e 2 / n 2 b. We 
can find Zi (T) by equating this energy with the decrease 
in the free energy per electron ksT ln(no/n) due to the 
entropy increase (the entropy increase can be derived ac¬ 
cording to § 104 of Ref. l32l . where ks is the Boltzmann 
constant, n = Zi(T)N is the concentration of ionized 
electrons and no = 2/A 3 with A = a/27 tH 2 /vPksT as 
the DeBroglie wavelength of free electrons at tempera¬ 
ture T. At k = 20000, b = 0.29 A, m* = 1.8?n e where 
m e is the electron mass and N = 10 15 cm -3 (estimated 
from that the concentration of total donor electrons is 
around 10 18 cm -3 and each cluster contributes ~ 300 
donor electrons), we get Zi{T) > Z* at T > 8 K with 
Z* = 100 which is a reasonable estimate. This shows 
that the outer electrons are completely ionized at tem¬ 
peratures that are not too low. For the inner core elec¬ 
trons, the dielectric response is nonlinear and the attrac¬ 
tive potential is stronger. So the ionization energy is 
higher ~ A(Z*e/A tt) 3 /5Pq R 5 for electrons at r ~ R. At 
R = 4a, it is found that only at T > 450 K can the in¬ 
ner electrons be ionized by a considerable quantity (the 
1/r 15 / 2 tail is completely stripped then). So the inner 
electrons are robust against the thermal ionization. 

Experimentally, charged clusters can be created con- 
trollably on the surface of LAO/STO structure when the 
LAO layer is of subcritical thickness < 3 unit cel l 47 * 48 * . A 
conducting atomic force microscope (AFM) tip is placed 
in contact with the top LaAlOa (LAO) surface and biased 
at certain voltage with respect to the interface, which is 
held at electric ground. When the voltage is positive, a 
locally metallic interface is produced between LAO and 
STO where some positive charges are accumulated in the 
shape of a disc. The same writing process can also create 
a periodic array of charged discs. 


Let us first concentrate on a disc of positive charge 
created in this manner on the STO surface. Close to the 
surface and in the bulk STO, one should apply the plane 
solution given by Eq. (14). When the distance r from 
the disc center is large, i.e., r R, the disc behaves 
like a charged sphere. Our results for a sphere are still 
qualitatively correct in this case. 


In a periodic array of highly charged discs with period 
2 L (see Fig. 3a in Ref. ITfl) . the linear concentration of 
free electrons responsible for the conductance at a very 
low temperature is of the order of n(L)L 2 , where n(r) 
is the electron density around a spherical donor cluster 
given by Sec. [Vi} When the overlapping parts between 
neighboring discs belong to the outer electron shells, the 
corresponding density at r = L is that of a Thomas- 
Ferrni atom with charge Z*. In this situation, the over¬ 
lapping external atmosphere forms conductive “bridges” 


between discs at low temperature. When T increases, 
however, the outer electrons are ionized and the bridges 
are gone. These free electrons spread out over the bulk 
STO. At T < 30 K, electrons are scattered mainly by 
the Coulomb potential of donors and the correspond¬ 
ing mobility decreases with a decreased electron velocity. 
For the electrons ionized into the vast region of the bulk 
STO, they are no longer degenerate, so their velocity be¬ 
comes much smaller at relatively low temperature. This 
results in a much smaller mobility of the ionized electrons 
than those bound along the chain. Their contribution to 
the conductivity is thus negligible. The system becomes 
more resistive due to the ionization and one can observe 
a sharp decrease of the conductivity along the chain. 


IX. CONCLUSION 

In this paper we have studied the potential and elec¬ 
tron density depth profiles in surface accumulation layers 
in crystals with a large and nonlinear dielectric response 
such as SrTiOa (STO) in the cases of planar, spherical 
and cylindrical geometries. We use the Landau-Ginzburg 
free energy expansion for the dielectric response of STO 
and adopt the Thomas-Fermi approximation for its elec¬ 
tron gas. 

For the planar geometry we predict an electron den¬ 
sity profile n(x ) oc (x + d) -12//7 , where d oc D Q 7 ^ 5 with 
Do as the surface induction and x as the distance from 
the interface. Here we skipped comparison of this result 
with the experimental data which can be found in Ref. 
EDI The data generally show a reasonable agreement with 
our predictions. We extend our studies from a single ac¬ 
cumulation layer to overlapping ones, and also investigate 
the electron “spill-out” from a heavily n-type doped STO 
(NSTO) into a moderately n-type doped STO. 

In the second part of this paper we study the collapse 
of electrons onto spherical and cylindrical donor clusters. 
Such “ fall-to-the-center” originates from the very fast 
decrease of the electron potential energy near the cluster 
which is oc (—1/r 5 ) in the spherical case and in turn is 
a result of the strongly nonlinear dielectric response of 
STO. This leads to a very unusual two-scale shape of the 
electron density around the cluster. We show how one 
can verify such a shape experimentally. 
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